10 



Journal of Modern Optics PQE'l 



Journal of Modern Optics 

Vol. 00, No. 00, 00 Month 200x, 1-8 



Spectroscopic Characterization of Gapped Graphene in the Presence of 

Circularly Polarized Light 

Godfrey Gumbs 1 , O. Roslyak 1 , Danhong Huang 3 , Antonios Balassis 2 

1 Department of Physics and Astronomy, Hunter College at the City University of New York, 695 Park 

Avenue New York, NY 10065, USA. 
2 Physics Department, Fordham University, 441 East Fordham Road, Bronx, NY 10458, USA. 
3 Air Force Research Laboratory, Space Vehicles Directorate, 
Kirtland Air Force Base, NM 87117, USA. 

(Received 00 Month 200x; final version received 00 Month 200x) 

We present a description of the energy loss of a charged particle moving parallel to a graphene layer and 
graphene double layers. Specifically, we compare the stopping power of the plasma oscillations for these two 
configurations in the absence as well as the presence of circularly polarized light whose frequency and intensity 
can be varied to yield an energy gap of several hundred meV between the valence and conduction bands. The 
dressed states of the Dirac electrons by the photons yield collective plasma excitations whose characteristics are 
qualitatively and quantitatively different from those produced by Dirac fermions in gapless graphene, due in 
part to the finite effective mass of the dressed electrons. For example, the range of wave numbers for undamped 
self-sustaining plasmons is increased as the gap is increased, thereby increasing the stopping power of graphene 
for some range of charged particle velocity when graphene is radiated by circularly polarized light. 
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1. Introduction 

Graphene [l[ has now been investigated intensively both theoretically and experimentally to 
gain a clearer understanding of its properties to see how it can be employed in nanoelectronic 
device applications. Several review articles have drawn together many aspects of graphene with 
regard to the role which many-body effects have on their electron transport and optical response 
properties 0-0] • However, we still need to explore in greater detail the energy loss when a beam of 
electrons travels near a single layer or double layers of graphene. Generally, inelastic interactions 
include phonon excitations, inter and intra-band transitions, plasmon excitations, inner shell 
ionizations, and Cerenkov radiation. This is a non-trivial challenge which we investigate by 
considering a model which accounts for only plasma losses and those due to intra-band and 
inter-band transitions. The beam of incident electrons is assumed to be focused in a narrow angle 
with a known, small range of kinetic energies. Some initial work has been reported in Refs. [5K3] 
for gapless graphene. In this manuscript, our interest is to investigate how the stopping power 
is affected when a gap opens up between the valence and conduction bands, as may occur, for 
instance, when graphene is radiated by a circularly polarized electric field (CPEF). 

Electron energy loss spectroscopy (EELS) has been considered for many systems dating back 
to Ritchie's classic paper 0] for a slab of dielectric material in the local limit and subsequently 
generalized to a non-local dielectric function by Gumbs and Horing (jj. The nonlocality was 
included in Ref. 0] through the screening of the electron-electron interaction produced by single- 
particle excitations and plasmon modes. Fessatidis, et al. [3] recently adopted the formalism of 



Horing, et al.jlOl] for a 2D electron gas and Gumbs and Balassis U| for a nanotube to the case 



of gapless graphene with the aid of the polarization function calculated by Wunsch, et al. 
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for conventional Dirac electrons. 

There have been several recent papers which considered the effects which a CPEF [l3l . 14]. 
spin-orbit interaction (SOI) [TH] in suspended graphene or the sublattice symmetry breaking 
(SSB) by an underlying polar substra te |16l -ll8l] in epitaxial graphene may have on the energy 
band structure and plasma excitations [14] of a graphene sheet. Under these conditions, there is a 
gap between the valence and conduction bands as well as between the intra-band and inter-band 
electron-hole continuum of the otherwise semimetal Dirac system. 1^, lil, 20] Additionally, the 
interplay between the single-particle excitations in the long wavelength limit results in dielectric 
screening of the Coulomb interaction which produces an undamped plasmon mode which appears 
in the gap separating the two types of electron-hole modes forming a continuum. By this, we 
mean that a self-sustained collective plasma mode is not supported by exciting either the intra- 
band or inter-band single-particle modes only. Furthermore, although the Dirac electrons near 
the K-points now acquire a non-zero effective mass with a CPEF, one still cannot produce a long 
wavelength plasmon mode as it is possible in the two-dimensional electron gas (2DEG) by intra- 
band excitations only. These properties of the plasma modes give rise to noticeable differences in 
the behavior of the stopping power of gapped graphene compared with conventional graphene. 

A dynamic energy gap produced by the CPEF may be comparable with that induced by a polar 
substrate. However, unlike the latter, the gap produced by a CPEF is highly tunable. Addition- 
ally, those two gap-opening mechanisms may take effect simultaneously. Whenever this occurs, 
the circular polarization of the electromagnetic field may no longer be required provided that it 
is in resonance with the SSB gap. Electron energy loss spectroscopy (EELS) may be employed 
to ascertain the plasmon frequencies in single and double layer graphene. This setup comprises 
an electron/photon mixed pump-probe, with the electron beam being the probe effected by the 
polarized light whose role is to open the gap. The Raman shift of the scattered electrons provides 
both particle- hole and plasmon excitation frequencies, which are usually characterized by their 
spectral weight, a quantity that depends on the transferred energy frw and momentum hq. 

The outline of the remainder of this paper is as follows. In Sec. [2j we present the model Hamil- 
tonian for Dirac electrons in the presence of a photon field. We present the energy eigenvalues 
and eigenfunctions of the dressed Dirac electrons in the presence of circularly polarized light. 
We present the RPA result for the polarization matrix for a double layer in Sec. [3] which we then 
employ to obtain the energy loss for a charged particle moving parallel to the graphene sheets. 
Our numerical results for the energy loss are presented in Sec. HI 



2. Dressed Dirac electrons and Floquet sub-bands. 

Let us consider several graphene layers separated by distance d from each other. The layers may 
be epitaxially grown on a carbon or silicon based polar substrate. In the latter case the system it 
to be embedded into a micro-cavity (MC) for effective electron dressing. The central MC mode is 
characterized by associated energy NqHujo, with the mode being tuned with the substrate induced 
gap. Here, iVo is the average number of photons in the mode, and depends on specifics of the 
excitation method. For example, if the mode is thermally excited then Nq ~ hujo/ksT for small 
thermal energy k^T . Another method of excitation is to use an external circularly polarized laser 
beam making Nq a function of the external pump intensity. This dressing method is suitable for 
gapless free standing graphene. In our formalism, specific details of the method of excitation do 
not play a crucial role and will be neglected. Additionally, we assume that all graphene layers are 
in the node of the optical filed and thus we neglect retardation effects. This is valid if d <C 2tt/ujq. 

In order to proceed further, we must specify the model Hamiltonian at the two inequivalent 
K and K' points in the graphene energy dispersion, i.e., 
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W 



a+a + cr_a 



Hvf ((o+ + cr_) k x ± i (cj_ — <r+) fc^ 



(1) 
(2) 
(3) 



The operators involved in the above Hamiltonian act in the joint electron-photon space via the 
following relations: 



a\ it, N ) 



ff+ = lt><4.| 

a + |t,iVo} = 
<r-\i,N )=0 
W \it,N -l) 



a- = \i)(t\; 

a + \i,N ) = \t,N ) 
(T_|t,iVo) = |^JVo> 
at | |t, No) = VWTT\ |t, No + 1) 



with |t|) being the Dirac pseudo-spin basis. In this notation, the Jaynes-Cummings (JC) 
Hamiltonian ([2]) takes into account the circularly polarized field of frequency ujq and ampli- 
tude sj AirNohwo/V. Also, V is the characteristic volume of the mode. The Dirac part of the 
Hamiltonian is given by Eq. ([3]), with ± corresponding to the K,K' points, respectively. 

The energy of an electron rotational motion induced by the field is denoted by Wo oc \/Nq. 
Provided it is assumed to be much less than the energy of the optical filed itself, i.e., Wo/fiwo <C 1, 
the Hamiltonian may be further simplified [HI] to: H = INqHujo + (E g /2) cr 3 + T-Ld where E g = 
hwo is the energy gap between valence and conduction bands. The gap provides 



dressed Dirac electrons with an effective mass m* = E g j (2h 2 v F ^ . Corresponding eigenvalues and 
eigenfunctions summarized in Table lATl Fl . For the reader's convenience, dressed and free Dirac 
electrons are compared. It is a simple matter to show that the latter are just the off-resonance 
limit of the dressed states. 
There is an alternative approach to the problem of the gap opening by light based on classical 



description of light [22B25J]. In their approach the quasi-energies of the optically induced sub- 
bands were obtained by exploiting the Floquet theorem (See Eq.(5) in Ref.[24]). When the laser 
field of intensity ~ Aq changes its polarization from being linear (0 = 0) to circular polarized 
((f) = pi/2) the degeneracy of the subbands (±hvpk + NqMY) is lifted at the crossing points. 
The dynamic gap occurs at the Dirac point (k = 0). It is proportional to the intensity of the 
optical field: E g = 2 sui((f))(eAoVF) 2 /Kl. The dynamic gap also occurs at the cross-points away 
from the Dirac point (k = frw/2vf). It is proportional to the amplitude of the optical field: 
Eg = eAoVF^l — cos(<p) sin(2a), where tan(a) = k y /k x . 

A few words must be said about the effective chemical potential for the Floquet states, which 
separates the occupied (hole-like) from the unoccupied (electron-like) states in quasi-equilibrium. 
Its role is played by the "mean" energy (See Eq.(13) in Ref.Q]), which depends on the frequency 
of the light. It will be used in the next section to define the distribution of the Floquet states 
(dressed Dirac electrons). 

In conventional semiconductors the optically induced side-bands are responsible for effects 
such as photon-assisted tunnelingj26l. \2l\ , generation of excitons[28, 29], Stark[2§] and dynamical 



^^Note that the dressing does not mix the vales. Only K point is rpresented in the table since the relevant parameters 
(transition energies and the structure factor) are valey independent. 
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Frantz-Keldish phenomena[30l. 31]. Ingraphene, the side bands formation was analyzed in con- 
nection with photo- voltaic Hall effect 22], dc transportflU and pump-probe optical response [23]. 
The next section focuses on yet another effect - electron energy loss via inelastic scattering by 
plasmons. 



3. Energy Loss Formalism 

Let us consider a charged particle with charge Ze moving parallel to the graphene surface with 
velocity v at a distance zq from the top layer of a pair of graphene sheets with separation d. 
Then, the rate of loss of energy of this charged particle due to the frictional force it experiences 
by passing over the interacting electron gas is given by 



— = -2*^ J dz J —to.v)—— shnr («,,*;g,-q-v) (4) 

where e _1 (z, z'\ q, ui) is the inverse dielectric function. For a double layer, we may follow the 
procedure presented in Ref. (2l| to obtain the inverse dielectric function as a single layer may 
be described by the same formalism by setting the inter-layer distance to infinity. Choosing the 
origin z = at one of the layers, we obtain the inverse dielectric function as [2ll] 



e-\z,z';q,u) = 5{z-z')+ £ 8 (z' - (j - l)d) v c (jd, /d)%,(g, u) (5) 

jj'=o,i 

where 5(z — z') is the Dirac delta function and polarization matrix is 



u u u 12 \ i /nff (i-. c (g)ng) , c (g) e-^nffng , 



n 21 n 22 ; e ( 9>w ) ^ Ue(g) e-^nffng n£> (1 - « c ( 9 )n<? 

It is expressed in terms of the generalized dielectric function 



e(q,u;) = l-v c ( q ){[u^+U^ 



+ 



i- e -^)t, e (g)ng)ng)]} , 



(7) 



where v c (q) = 2ne 2 /e sq with e s = Aneoeb, e& being the average background dielectric constant. 
Noninteracting polarization of the given graphene layer assumes Lindhard form 



tt(0)(„,,\ - 1 V" /o(^k,A) ~ /o(^k+q,A') r AA' /on 

11 (9 '" )ll(22) "A,. ^ „ ^-^.v+^a k ' k+q ' (8) 



k, A,A'=±1 



Here, we have E^x = Xhvpk for free Dirac electrons, but \ = hvp ^Nqojq + X^J (e g /2) 2 + fc 2 | 

for the dressed ones. Also, the corresponding structure factors J^ k+ and particle distributions 
/o(-^k a) are given in Table [AT] A is the normalization area in the reciprocal space chosen to 
assure (1/A) « J <i 2 k. The upper limit of the integration is given by the ultraviolet cut-off as 
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in Ref . [121 ] . Straightforward substitution of Eq. ([5]) into (jlj) yield the energy loss for graphene 
double-layer: 

dW „ {Zef f d 2 q_ (q-v) . . 

-— — = F • vii = 2-7T — — — / - — To -v c (q)x 9) 

dt 11 4vre i (2vr) 2 g w v ; 

x9m ^e- 2qZo U u {q, -q • v) + e- 2q{zo+d) U 22 (q, -q • v) + 2e- q{2za+d) U l2 (q, -q • v)} 

The single layer stopping power may be obtained from Eqs. (|6l7l9p by putting 11^ (q, w) = 0. 

The contributions to the polarization matrix and consequently to the rate of loss of energy in 
the above Eq. may be identified as coming from two terms 0, 10, [IH. One is determined by 
the plasmon excitations, i.e., when both the real and imaginary parts of e(q, oj) are simultaneously 
equal to zero. The other comes from single-particle excitations, i.e., whenever 9m e(q,ui) is non- 
zero. Their relative contribution as a function of the energy gap induced by the dressing is the 
subject of the next section. 



4. Numerical results and Discussion 

In Fig. [U we plot the rate of loss of energy of a charged particle moving with speed v at distance 
zq = 2kp 1 from a single and double graphene layer configuration. For the pair of layers, we 
assume that the energy gap is the same for both layers whose separation we chose as d = kp 1 . 
Both layers are equally doped to yield the chemical potential energy fi, which is related to 
the Fermi level as \J [i 2 + (E g /2) 2 = Uvf^f- To evaluate the noninteracting polarization (|8|) all 
frequencies in the polarization function are analytically continued into the complex plane and 
acquire a small positive imaginary part so that uj — y co + i'y. In that case the polarization assumes 
the form: 

u( q , u + vy) = in -— ^nm Up \\ 6 {UJ ~ ^ + n(9 ' w) (10) 

{d/du)Qmm > (q, oj p ) \ u==Up 

where we have omitted the layer indices for brevity. The plasmon resonances, denoted as lj p , 
correspond to the poles of the inverse dielectric function which are given by the solutions of 
e(q,ojp — i*y) = 0. For infinitesimally small 7 = + an analytical solution is readily available 
(See [14| and references therein). We chose j/fj, = 10~ 4 to account for possible lifetime broad- 
ening. The noninteracting polarization was calculated via adaptive local Monte-Carlo numerical 
integration. The energy loss simulation is plotted in Figs. [T]for various values of E g . 

The simulations shows that the plasmon contribution to the energy loss is decreased with 
growing E g , while the particle- hole contribution remains relatively unaffected by the gap as 
illustrated by Fig. [TJ This may be attributed to diminishing of the plasmon frequency u) p (q) ~ 
\J q(l — (E g /2fi) 2 ) as shown in Fig. [2] (a). For the double layer configuration, a small spike at 
low velocity in the plasmon contribution to the energy loss is due to formation of the undamped 
acoustic (low frequency) and optical (high frequency) plasmon branches. The plasmon modes 
whose frequency is within the regions 



PH> = Huj > Jh 2 v 2 F q 2 + E 2 (11) 



PH< = < huj < fi + ^h 2 v 2 (q + k F ) 2 + {Eg/2) 2 (12) 
are Landau damped and given by $tee(q,ujp — 17) = 0, and also depicted in Fig. [2j The damped 
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Figure 1. (Color online) Stopping power in units of (Ze/2ir) 2 eJ 1 as a function of charged particle velocity in units of 
vf- Panels (a) correspond to double layers, panels (b) to single layer configuration. Rows (1), (2) and (3) correspond to 
Eg/fi = {0.0,1.0,1.5}, respectively, and are depicted by the red, green and blue curves, correspondingly. The solid black 
curves are the particle-hole contributions, and the dashed curves represent plasmon contributions. Both (a. 4) and (b.4) 
compare particle hole mode contributions for the three values of the chosen energy gaps, and (a. 5) and (b.5) compare the 
corresponding plasmon contributions. 




Figure 2. (Color online) Undamped (two higher frequency branches) and damped (two lower in frequency branches) plasmon 
dispersion relations for gapped graphene. The curves in (a) correspond to a single layer, the curves in (b) to a double layer 
configuration. The red, green and blue curves show the plasmon dispersion for E g //i = {0.0, 1.0, 1.5}, respectively. 
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plasmons contribute as particle-hole excitations. 

The plasmon contribution to the energy loss may be qualitatively explained by spanning the 
energy loss sector < ui < qv through the plasmon dispersion if Fig. [TJ In the low velocity regime, 
given by PH<(u; — > qv), one has a contribution from the intraband particle- hole excitations 
which are relatively weak due to the suppressed backscattering. Once v is increased so that 
the plasmons start to contribute one has rapid increase in the lost energy. Once the sector fully 
engulfs the plasmon mode (modes) their contribution reaches its maximum and becomes less than 
the particle- hole mode contribution when v/vp ~ 8. After that the increase in the energy loss 
is governed by the strong interband particle-hole excitations of PH> . Both types of the particle- 
hole excitations are given by the second term in Eq. (jlOp . For gapless graphene the undamped 
plasmon may live only in the triangular region where the distractive interference between the two 
types of the particle-hole excitation reduce the imaginary part of the polarization to zero. The 
gap provides an additional space for the plasmons to run to the higher values of the wave vector. 
In the double layer configuration both optical and acoustic plasmon modes contribute, but as 
v is increased, only the optical mode contributes. Since the spectral weig ht0 of the acoustical 
mode is large compared with the optical branch, we see a spike once its contribution reaches the 
maximum. 

The effect of the energy gap on the plasmon contribution to the energy loss is summarized in 
Figs. |2](a.5) and[2]Jh.5). These results show intricate interplay between the plasmon and particle- 
hole contributions to the energy loss, governed by the window in the particle-hole continuum is 
opened up in the presence of CPEF. On one hand the gap reduces the velocity threshold of the 
collective plasma excitations. The shorter wave length plasmon frequency become available in 
the window, and may be accommodated by smaller energy loss sector. Furthermore, the plasmon 
branch decreases in frequency as E g is increased, as seen in Fig. [21 thereby making it more likely 
to excite plasmons even when the charged particle speed is low, and making gapped graphene 
to have a higher stopping power than conventional gapless graphene. On the other hand, the 
plasmon contribution reaches its maximum for smaller velocities as and the role these collective 
and particle-hole plasma excitations play in the energy loss are reversed in the high-velocity 
limit. In this regime where the particle-hole excitations should dominate the window makes the 
stopping power of gapped graphene being less than when there is no gap. 
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the factor in front of the delta function in Eq. (I X 1 ) . 
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Tabic Al. Free vs dressed Dirac electron states. 



Free Dirac 



Dressed Dirac 



E k = N oj ± ^(£ g /2) 2 + fc 2 



Energy 



E k = ±k 



f N ) 



Wave 
Function 



1 1 



1 e l 



It) 
ID 



1 aAfo,fce 



2i- 



k,k+q - (l+o 2 W0 , J(l+^ D , fc+ ,) ^ + «7V0,fc«JV0,fe+g- 
±2aAr 0j A:aArO,A;+g cos &k,k+q 



Structure 
Factor 



•^k.k+q = f(l ±cos6» fcjfc+g ) 



P_ exp-(g k -ri)/fc B :r~ 
k l+exp-(£ k -^)/fc B T 



f/3 _ cxp -{E±-n-Nab> )/k B r r~ 
Jk 1+cxp -{E k -fi-N a u> )/k B T 



Distribution 



Parameters 



9 fc = tan- 1 (j* 



awo.fc 



It) 



COS t^fc+q 



;ll) 



l + |^o) 

\-lNo) 



^k 2 +(E g /2) 2 +E g /2 



COS 



sin 



ItJVb) 
l4-,JVb> 



fc+g cos ( 



A /fc 2 +g 2 +2fcgcose fc ,, 



tan 2 (|) 



Va;g+4iyo 2 (Afo + l)/Afo+^o 
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